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Heat  transfer  through  the  gas  diffusion  layer  (GDL)  is  a  key  process  in  the  design  and  operation  of  a 
proton  exchange  membrane  (PEM)  fuel  cell.  The  analysis  of  this  process  requires  determination  of  the 
effective  thermal  conductivity.  This  transport  property  differs  significantly  in  the  through-plane  and 
in-plane  directions  due  to  the  anisotropic  micro-structure  of  the  GDL. 

A  novel  test  bed  that  allows  separation  of  in-plane  effective  thermal  conductivity  and  thermal  contact 
resistance  in  GDLs  is  described  in  this  paper.  Measurements  are  performed  using  Toray  carbon  paper  TGP- 
H-120  samples  with  varying  polytetrafluoroethylene  (PTFE)  content  at  a  mean  temperature  of  65-70  °C. 
The  measurements  are  complemented  by  a  compact  analytical  model  that  achieves  good  agreement 
with  experimental  data.  The  in-plane  effective  thermal  conductivity  is  found  to  remain  approximately 
constant,  17.5  Wm-1  K-1,  over  a  wide  range  of  PTFE  content,  and  its  value  is  about  12  times  higher 
than  that  for  through-plane  conductivity. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  temperature  distribution  in  a  proton  exchange  membrane 
(PEM)  fuel  cell  is  non-uniform  due  to  the  electrochemical  reac¬ 
tion  and  associated  irreversibilities  [1-5].  Accurate  knowledge  of 
the  temperature  distribution  and  associated  heat  transfer  rates  is 
required  to  determine  various  transport  phenomena  such  as  water 
and  species  transport,  reaction  kinetics,  and  rate  of  phase  change. 
Thermal  transport  also  impacts  design,  efficiency,  reliability  and 
durability  of  the  system  [6-8]. 

A  key  thermo-physical  property  for  thermal  analysis  of  fuel  cells 
is  the  thermal  conductivity  of  the  membrane-electrode  assembly 
components,  particularly  the  gas  diffusion  layer  (GDL)  [9,10].  The 
fibrous  anisotropic  micro-structure  of  a  GDL  combined  with  the 
large  differences  between  the  thermal  conductivity  of  the  solid 
(carbon  fibers)  and  fluid  (air/water)  phases  make  it  challenging  to 
determine  the  dependence  of  the  effective  thermal  conductivity  on 
direction. 

The  majority  of  the  thermal  analyses  of  fuel  cells  have  relied 
on  a  simplified  model  representation  that  assumes  an  isotropic 
thermal  conductivity  [10-12]  that  is  determined  as  a  combination 
of  the  parallel  and  series  models  and/or  based  on  the  geometric 


*  Corresponding  author  at:  Department  Mechanical  Engineering,  and  Institute  for 
Integrated  Energy  Systems,  University  of  Victoria,  P.O.  Box  3055,  Victoria,  BC,  V8W 
3P6,  Canada.  Tel.:  +1  7787828587;  fax:  +1  2507216051. 

E-mail  address:  ehsans@uvic.ca  (E.  Sadeghi). 

0378-7753 /$  -  see  front  matter  ©  2010  Elsevier  B.V.  All  rights  reserved, 
doi:  10.101 6/j.jpowsour.2010.1 1.1 51 


mean  of  the  thermal  conductivities  of  solid  and  fluid  phases  [11,12]. 
Although  the  thermal  resistance  of  GDLs  for  the  in-plane  direction 
is  higher  compared  to  the  through-plane  direction,  heat  transfer  to 
the  bipolar  plate  (BPP)  occurs  in  both  directions  due  to  the  alternat¬ 
ing  nature  of  the  land  and  channel  areas  [13].  In  the  few  modeling 
studies  that  have  considered  anisotropy,  parametric  investigations 
have  shown  that  the  prescription  of  anisotropic  properties  has  a 
major  impact  on  current  density  distribution  and  on  the  relative 
importance  of  limiting  transport  processes  [  1 3-1 5  ].  The  determina¬ 
tion  of  the  in-plane  thermal  conductivity  is  therefore  an  important 
parameter  for  thermal  analysis  and  management  of  PEM  fuel  cells 
and  stacks. 

Theoretical  prediction  of  the  in-plane  thermal  conductivity  of 
GDLs  includes  the  work  of  Zamel  et  al.  [1 6]  who  developed  a  numer¬ 
ical  model  to  estimate  the  through-plane  and  in-plane  effective 
thermal  conductivities  of  a  dry  untreated  carbon  paper  GDL.  They 
studied  the  effects  of  porosity,  fiber  distribution  and  compres¬ 
sion  on  the  effective  thermal  conductivity  and  concluded  that  the 
impact  of  fiber  distribution  is  more  pronounced  for  the  through- 
plane  direction  than  the  in-plane  direction.  The  numerical  results 
indicate  that  porosity  is  an  essential  determinant  of  the  effective 
thermal  conductivity  of  a  GDL  but  not  compression.  Based  on  the 
results,  Zamel  et  al.  [16]  proposed  correlations  for  the  effective 
thermal  conductivity  of  a  dry  GDL  with  no  binder  or  hydrophobic 
treatment. 

GDLs  are  generally  treated  with  polytetrafluoroethylene  (PTFE) 
to  render  them  hydophobic  and  enhance  liquid  water  transport 
[17].  The  effect  of  PTFE  treatment  on  the  thermal  transport  param- 
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Nomenclature 

A 

cross-sectional  area  of  fluxmeters,  heat  radiation 
area  (m2) 

Af 

cross-sectional  area  of  a  fiber  (m2) 

An 

cross-sectional  area  of  GDL  (m2) 

a 

radius  of  contact  area  between  fibers  (m) 

b 

radius  of  the  area  covered  by  PTFE  at  contact  points 

(m) 

e 

emissivity,  Eq.  (1) 

Fij 

view  factor,  Eq.  (1) 

k 

thermal  conductivity  (W  m-1  K-1 ) 

^cell 

in-plane  effective  thermal  conductivity  of  each  cell, 
Fig.  3  (Wm-1  K-1) 

^eff.in 

in-plane  effective  thermal  conductivity  of  GDL 
(Wm-’K4) 

^eff.in.o 

in-plane  effective  thermal  conductivity  of  GDL,  no 
PTFE(Wm_1  K-1) 

^parallel 

effective  thermal  conductivity  based  on  the  Parallel 
model  (Wm-1  K-1) 

kpiFE 

thermal  conductivity  of  PTFE  (W  m-1  K-1 ) 

ks 

thermal  conductivity  of  carbon  fiber  (W  m-1  K-1 ) 

L 

length  of  GDL  sample,  distance  between  sample 
holders,  Figs.  1  and  3  (m) 

m 

number  of  cells 

N 

total  number  of  GDL  samples  in  the  experiment 

n 

number  of  fibers  in  each  cell 

Q 

heat  transfer  rate  (W) 

Qij 

Radiation  heat  transfer  between  bodies  i  and  j  (W) 

Q 

heat  flux  (Wm-2) 

RceW 

cell  thermal  resistance  (K  W-1 ) 

^end 

thermal  resistance  at  sample  ends  (K  W-1 ) 

Rfl 

thermal  resistance  of  fluxmeter  (K  W-1 ) 

Rgdl 

GDL  thermal  resistance  (K  W-1 ) 

Rgr 

groove  thermal  resistance  (I<  W-1 ) 

A  unc 

total  thermal  resistance  at  each  cell  interface 
(KW-1) 

Aayer 

summation  of  GDL  and  groove  thermal  resistance 
(KW-1) 

RpiFE 

thermal  resistance  of  PTFE  in  contact  regions  of 
fibers  (KW-1) 

KsH 

sample  holder  thermal  resistance  (KW-1 ) 

^tot 

total  thermal  resistance  (KW-1 ) 

Rw 

thermal  resistance  of  wooden  block  (KW-1 ) 

r 

fibers’  mean  radius  (m) 

T 

temperature  (K) 

TCR 

thermal  contact  resistance  between  fibers  (KW-1 ) 

t 

sample  thickness  (m) 

V 

average  volume  of  a  fiber  (m3) 

Vs 

total  volume  of  fibers  (m3) 

Vtot 

total  volume  of  GDL  sample  (m3) 

W 

width  (m) 

z 

heat  flow  direction  (m) 

Greek  symbols 

a 

radius  of  contact  area  between  fibers  over  fiber 
radius,  ar-2 

p 

radius  of  area  covered  by  PTFE  over  fiber  radius,  hr-2 

X 

weight  fraction  of  PTFE  content 

G 

Stefan-Boltzmann  constant 

so 

GDL  porosity  before  PTFE  treatment 

0 

angle  between  fiber  and  heat  flux  directions 

An 

maximum  fiber  angle 

0s  fiber  volume  fraction 

l  fiber  average  length  (m) 

Subscripts 

1  right  side  of  the  experiment  setup 

2  left  side  of  the  experiment  setup 

i  body  i 

j  body j 

low  lower  fluxmeter 

r  radiation 

up  upper  fluxmeter 


eters  has  been  investigated  by  Khandelwal  and  Mench  [18].  Their 
measurements  for  SIGRACET®  GDLs  showed  that  for  20%  PTFE  con¬ 
tent,  the  through-plane  thermal  conductivity  was  reduced  by  54% 
compared  to  untreated  GDL  samples;  the  thermal  contact  resis¬ 
tance  on  the  other  hand  was  not  significantly  affected  by  the 
variation  of  PTFE  content.  These  results  are  opposite  to  expected 
trends  from  physical  considerations;  additional  PTFE  is  for  instance 
expected  to  displace  lower  conductivity  air,  and  hence  results  in 
higher  effective  thermal  conductivity.  Karimi  et  al.  [  1 9]  investigated 
the  effect  of  PTFE  coating  on  the  through-plane  conductivity  and 
thermal  contact  resistance  of  SpectraCarb  GDLs,  and  reported  that 
the  effective  thermal  conductivity  of  PTFE-treated  GDLs  increased 
slightly  at  low  compression  loads  and  decreased  slightly  at  higher 
loads.  For  low  compression  loads,  they  also  reported  significantly 
higher  thermal  contact  resistance  values  for  PTFE-treated  samples 
compared  to  untreated  ones.  This  difference  decreased  with  an 
increase  in  applied  pressure. 

Several  experimental  approaches  have  been  proposed  to  mea¬ 
sure  electrical  conductivity  of  GDLs.  Ismail  et  al.  [20]  measured  the 
in-plane  and  through-plane  electrical  conductivities  of  SIGRACET® 
GDLs,  and  found  the  in-plane  electrical  conductivity  remains 
approximately  constant  with  an  increase  in  PTFE  content.  They 
also  reported  that  the  through-plane  electrical  contact  resistance 
increases  with  PTFE  content. 

The  micron  scale  of  the  fibers  combined  with  the  brittle  nature 
of  the  GDL  structure  make  it  challenging  to  measure  the  in-plane 
thermal  conductivity  of  such  random  micro-structures,  and  to  the 
authors’  knowledge,  no  experimental  data  has  been  reported  in  the 
open  literature. 

This  paper  presents  a  combined  experimental  and  theoretical 
investigation  focusing  on  the  determination  of  the  in-plane  ther¬ 
mal  conductivity  of  PTFE-coated  GDLs.  Building  on  our  previous 
study  that  dealt  with  the  through-plane  conductivity  and  contact 
resistance  [21],  the  existing  test  bed  was  modified  and  an  exper¬ 
imental  technique  was  developed  that  enables  the  measurement 
of  in-plane  thermal  conductivity  of  fibrous  porous  media  and  thin 
films.  Toray  carbon  papers  TGP-H-1 20  with  different  PTFE  contents 
are  used  in  the  experiments.  The  in-plane  effective  thermal  conduc¬ 
tivity  and  contact  resistance  are  deduced  from  the  total  thermal 
resistance  measurements  by  performing  a  series  of  experiments 
with  GDL  samples  of  different  lengths  but  similar  micro-structures. 
Furthermore,  a  compact  analytical  model  is  proposed  to  predict  the 
in-plane  thermal  conductivity  of  GDLs  as  a  function  of  porosity  and 
PTFE  content. 

2.  Experimental  study 

The  experimental  apparatus  and  a  schematic  of  the  test  setup 
for  the  in-plane  thermal  conductivity  measurement  are  shown  in 
Fig.  1.  The  test  chamber  consists  of  a  stainless  steel  base  plate  and 
a  bell  jar  enclosing  the  test  column.  The  test  column  consists  of, 
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Fig.  1.  (a)  Experimental  apparatus  used  for  in-plane  thermal  conductivity  measurement  and  (b)  schematic  view  of  the  test  column. 


from  top  to  bottom:  the  loading  mechanism;  the  heater  block;  the 
upper  fluxmeter;  the  upper  wooden  block;  the  sample  holders;  the 
samples;  the  lower  fluxmeter;  the  lower  wooden  block;  the  heat 
sink  (cold  plate);  and  the  polymethyl-methacrylate  (PMMA)  layer. 

The  heater  block  was  made  of  aluminum  in  which  a  pencil-type 
electrical  heater  was  installed.  The  designed  cold  plate  consisted  of 
a  hollow  copper  cylinder,  1 .9  cm  high  and  15  cm  diameter.  Cooling 
was  accomplished  using  a  closed  loop  water-glycol  bath  in  which 
the  coolant  temperature  can  be  set.  The  cold  plate  was  connected 
to  the  chiller  unit  which  adjusts  the  cold  water  temperature.  A  load 
was  applied  on  the  upper  wooden  block  and  fluxmeter  to  improve 
the  contact  between  the  sample  holders  and  the  fluxmeters. 

The  fluxmeters  were  made  of  standard  electrolyte  iron.  To  mea¬ 
sure  temperatures  along  the  fluxmeters,  six  T-type  thermocouples 
were  attached  to  each  fluxmeter  at  specific  locations  shown  in 
Fig.  1(b).  The  thermal  conductivity  of  the  fluxmeters  was  known 
and  used  to  measure  the  heat  flow  rate.  The  sample  holders 
were  made  of  aluminum  and  have  two  grooves  with  the  width  of 
1.15  mm.  Two  T-type  thermocouples  were  attached  to  each  sample 
holder  near  the  grooves  to  measure  the  temperature. 

2.1.  Sample  preparation 

Toray  carbon  papers  TGP-H-120  with  the  base  porosity  (poros¬ 
ity  of  the  untreated  GDL)  of  78%  and  a  wide  range  of  PTFE  content, 
from  5%  to  30%,  were  used.  The  thickness  of  GDL  samples  was  mea¬ 
sured  using  a  Mitutoyo  digital  micrometer  with  the  accuracy  of 
0.001  mm.  The  measurements  were  performed  10  times  for  each 
sample  at  different  locations,  and  the  average  values  are  reported 
in  Table  1 .  Rectangular  test  samples  were  cut  with  a  width  of  35  mm 
and  different  lengths. 


Table  1 

Thickness  of  examined  Toray  carbon  papers  TGP-H-120. 

30 

0.354 


2.2.  Test  procedure 

The  experiments  were  performed  under  a  vacuum  to  ensure 
negligible  convection  heat  transfer.  Depending  on  the  thickness, 
a  number  of  similar  GDL  samples,  e.g.,  three  sheets  for  TGP-H-120, 
were  stacked  together  and  inserted  in  each  groove  of  the  sample 
holders  as  shown  in  Fig.  1.  The  use  of  several  layers  of  GDLs  miti¬ 
gated  some  of  the  experimental  challenges  and  uncertainties:  use  of 
a  single  GDL  layer  between  the  sample  holders  leads  to  an  excessive 
temperature  drop  across  the  sample  holders,  which  in  turn  tremen¬ 
dously  increases  the  temperatures  along  the  upper  fluxmeter,  even 
for  a  small  heat  flux,  and  may  burn  the  thermocouples.  This  is  due 
to  micron  size  cross-sectional  area  of  the  GDL.  After  investigating 
and  trying  different  methods,  the  present  sample  holders  featur¬ 
ing  multiple  grooves  that  can  hold  several  GDLs  was  devised  to 
overcome  these  challenges. 

To  reduce  the  contact  resistance  between  the  groove  walls  and 
the  samples,  a  thin  layer  of  thermal  paste  was  applied  inside  each 
groove.  To  improve  the  stability  of  the  sample  holders  and  provide  a 
good  contact  with  the  fluxmeters,  compressive  loads  were  applied 
to  the  upper  wooden  block  and  the  upper  fluxmeter;  this  is  solely 
to  keep  the  test  column  together.  Thermal  paste  was  also  used  to 
reduce  the  thermal  contact  resistance  at  the  interfaces  between  the 
sample  holders  and  the  fluxmeters. 

Temperatures  were  monitored  continuously  and  recorded  when 
steady-state  conditions  were  achieved.  This  took  approximately 
7  h  for  each  experiment.  The  fairly  long  equilibration  time  is  due 
to  the  restricted  cross  sectional  area  through  which  heat  transfer 
takes  place.  The  temperature  gradient  between  the  hot  and  cold 
plates  results  in  one-dimensional  heat  conduction  from  the  top 
to  the  bottom  of  the  test  column.  The  thermal  resistance  network 
corresponding  to  the  experimental  setup  is  shown  in  Fig.  2(a).  Nat¬ 
ural  convection  heat  losses  are  negligible  in  the  vacuum  chamber. 
Radiation  heat  losses  from  the  fluxmeters  and  end  plates  can  be 
estimated  from  the  following  relationship  [22]. 

ofjf  -  If) 

^  ~  (1  -  e,/eA)  +  PP  +  (1  -  ej/ejAj) 


PTFE  (%) 
t  (mm) 


5 

0.374 


10 

0.376 


20 

0.362 


(i) 


3568 


E.  Sadeghi  et  al.  /  Journal  of  Power  Sources  196(2011)  3565-3571 


Fig.  2.  (a)  Complete  and  (b)  reduced  resistance  network  associated  with  the  experiments. 


where  a  =  5.67  x  10-8Wm-2  K-4  is  the  Stefan-Boltzmann  con¬ 
stant.  Qij  is  the  radiation  exchange  between  bodies  i  and  j,  and  e 
is  the  emissivity.  Also,  Fy  is  the  view  factor  defined  as  the  frac¬ 
tion  of  the  radiation  that  leaves  A*  and  is  intercepted  by  Aj.  To  find 
the  radiation  heat  losses,  Eq.  (1)  is  employed  which  provides  the 
maximum  radiative  heat  transfer  between  two  bodies.  The  inves¬ 
tigated  radiation  losses  from  the  fluxmeters  to  the  wooden  blocks 
and  the  chamber  wall  are  less  than  1%  of  the  total  heat  flow  pass¬ 
ing  the  fluxmeters.  Also,  small  temperature  difference  between  the 
fibers  in  GDLs  as  well  as  relatively  low  temperature  levels  (less  than 
370  K)  inside  the  medium  ensure  negligible  radiation  heat  transfer 
in  the  GDL.  Thus,  heat  transfer  is  only  due  to  conduction  and  can 
be  determined  using  Fourier’s  equation. 


where  dT/dz  is  the  temperature  gradient  along  the  test  column,  k 
is  the  thermal  conductivity  of  the  fluxmeters,  and  A  is  the  cross- 
sectional  area  of  the  fluxmeters.  Considering  negligible  heat  losses, 
the  resistance  network  shown  in  Fig.  2(a)  can  be  reduced  to  Fig.  2(b). 
The  total  thermal  resistance  between  two  sample  holders,  Rt0 1, 
includes  the  samples’  thermal  resistance  and  the  resistances  at 
the  sample  ends  (a  combination  of  the  thermal  contact  resistance 
between  the  grooves  and  the  samples  and  other  possible  resis¬ 
tances  caused  by  the  edges  of  the  grooves  at  each  end),  and  can 
be  expressed  as: 


thermal  conductivity. 

.  Ll  L2 

eff'ln  RcduA  RcduA 

(4) 

.  N(Li~L2 ) 

eff'm  (Rtotl  -Rtot2)An 

(5) 

where  l\  and  I2  are  the  sample  length,  the  distance  between  the 
two  sample  holders,  in  experiment  1  and  2,  and  Ain  is  the  in-plane 
cross-section  of  each  sample. 


2.3.  Uncertainty  analysis 

Considering  the  relationship  for  evaluating  the  in-plane  effec¬ 
tive  thermal  conductivity,  i.e.,  Eqs.  (3)  and  (5),  the  relevant 
parameters  in  the  analysis  can  be  expressed  as: 

/<eff,in=/(Q,AT,t,W,L)  (6) 

The  main  uncertainty  in  our  experiments  is  due  to  errors  in 
determining  the  heat  flux  through  the  sample  holders  which  leads 
to  a  maximum  error  of  3.7%.  The  maximum  uncertainties  for  the 
thermocouples  and  the  data  acquisition  readings  are  ±1  °C  which 
introduces  a  maximum  error  of  1.8%  between  two  sample  hold¬ 
ers.  Other  uncertainties  including  those  associated  with  the  width, 
thickness,  and  length  measurements  are  0.3%,  0.3%,  and  0.9%, 
respectively.  The  maximum  uncertainty  for  the  thermal  resistance 
measurements  can  be  calculated  from  [23]: 


n  ^layer  ^GDL  +  ^grl  +  ^gr2  RgDL  .  ^end  AT 

Rtot  =  ~N~  = - N - =  ~N~  +  ~N~  =  ~Q~"  (3) 

where  AT  is  the  temperature  difference  between  the  two  sample 
holders  and  N  is  the  total  number  of  GDL  layers  stacked  in  the 
grooves.  Agdl  and  Renc \  are  the  thermal  resistance  of  each  sample 
and  the  total  thermal  resistance  at  the  end  points  of  each  sample, 
respectively.  There  is  a  small  difference  between  the  heat  flow  val¬ 
ues  measured  for  the  upper  and  lower  fluxmeters  due  to  heat  losses 
to  the  lower  wooden  block  and  to  experimental  uncertainties;  heat 
losses  to  the  wooden  block  are  about  4%  of  the  heat  flow  passing 
the  upper  fluxmeter.  Therefore,  the  actual  heat  flux  which  passes 
through  the  GDL  samples  is  the  heat  flux  measured  at  the  lower 
fluxmeter.  To  ensure  accuracy,  this  heat  flow  rate  was  used  in  the 
analysis,  i.e.,  Eq.  (3).  To  find  the  in-plane  thermal  conductivity,  two 
sets  of  experiments  were  performed  with  different  sample  lengths. 
LJnder  the  same  experimental  conditions,  ftend  for  both  experiments 
was  assumed  to  be  equal.  Applying  Eq.  (3)  to  both  of  the  measure¬ 
ments  and  subtracting  them,  one  can  find  the  in-plane  effective 


(7) 


For  the  present  study,  the  maximum  uncertainty  is  estimated  to 
be  ±4.2%. 

3.  Analytical  study 

The  complex  micro-structure  and  associated  heat  transfer 
mechanism  of  fibrous  GDLs  make  it  difficult  to  develop  an  analytic 
model  for  the  effective  thermal  conductivity.  To  model  the  in-plane 
effective  thermal  conductivity,  a  random  micro-structure  divided 
into  m  equally-sized  cells  is  considered.  Each  cell  consists  of  n  fibers 
with  an  average  radius  of  r  and  an  average  length  of  l  which  are 
randomly  oriented  in  the  xy  plane  with  an  angle  0  to  the  in-plane 
heat  flow  direction  and  stacked  vertically  in  z  direction,  as  shown 
in  Fig.  3.  The  fiber  angle  0  can  vary  in  this  representation. 
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Fig.  3.  Constructed  micro-structure  for  in-plane  thermal  conductivity  modeling. 


Considering  that  the  primary  path  for  the  heat  conduction  is 
through  the  fibers  and  heat  transfer  between  the  fibers  in  a  cell 
is  negligible  due  to  large  contact  resistances,  a  parallel  equivalent 
circuit  model  can  be  used  to  determine  the  thermal  conductivity  in 
each  cell.  The  heat  conducted  through  the  ith  fiber  with  an  angle 
Oi  to  the  heat  flux  vector  is  cos  0 ,  therefore,  each  fiber  deviates 
from  the  parallel  model  by  cos  0[.  Thus,  the  effective  conductivity 
for  each  cell,  /<cen,  can  be  written  as: 


^cell  — 


y" 1  cos  o 


'^parallel  : 


E"=i(Mf  cose,) 

{nAf/<p s) 


=  cos  Ocpsks 


(8) 


where  cos  0  is  the  average  of  cos  Oi  in  a  cell.  ks  and  Af  is  are  the 
fiber  thermal  conductivity  and  cross-sectional  area,  respectively. 
0s  is  the  fiber  volume  fraction,  i.e.,  0S  =  1  -  £o,  where  £o  is  the  GDL 
porosity  before  PTFE  treatment.  We  assume  that  the  heat  flow  is 
transferred  from  cell  to  cell  through  the  junctions.  The  length  of 
each  cell  is  defined  as  the  average  conduction  path  in  a  cell,  ic os  6. 
Since  the  fibers  are  stacked  together  in  a  packed  micro-structure, 
it  is  assumed  that  each  fiber  of  two  neighboring  cells  is  in  contact 
with  two  fibers  from  the  top  and  the  bottom  and  carries  heat  from 
them  as  shown  in  Fig.  4.  The  contact  between  the  fibers  at  a  junction 
is  shown  in  Fig.  5. 

To  estimate  the  in-plane  effective  thermal  conductivity  of  the 
medium,  the  total  resistance  is  needed  which  can  be  found  using 
the  thermal  resistance  network  shown  in  Fig.  6. 

The  number  of  cells  as  well  as  the  number  of  fibers  in  each  cell 
is  required  to  evaluate  the  thermal  resistances.  These  values  can  be 


Fig.  5.  Contact  between  fibers  at  a  junction. 
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Rcell.l  \  Rjunc.l/  Rjunc,m-1  Rcell.m 

Fig.  6.  Thermal  resistance  network  of  the  GDL  for  in-plane  direction. 


found  through: 


L 

m  = 

t  cos  0 


(9) 


=  0sVW  =  Wt0scos6> 
mV  mV  Jtr 2 


where  V ,  W  and  t  are  the  average  volume  of  each  fiber,  the  width  and 
the  thickness  of  the  GDL  sample.  The  contact  resistance  between 
two  fibers  can  be  expressed  as  the  summation  of  constriction  and 
spreading  resistances  [24]. 


TCR  = 


2  ksa 


1 

2  ksra 


(11) 


where  a  is  the  radius  of  contact  area  between  fibers  and  a  is  a/r. 
A  portion  of  PTFE  in  GDLs  covers  the  contacting  fibers  providing 
an  additional  path  for  the  heat  flow  from  one  fiber  to  another.  The 
thermal  resistance  of  this  path  based  on  the  geometry  shown  in 
Fig.  5  can  be  expressed  as: 


1 

RpiFE 


RpiFE  = 


1  f  1  1  f  /<pxfe(27TX  dx) 

2j  dR  =  2  Ja  (r-y) 

1  (  udu 

^FTFEr  yja  1  _  0  _U2 


(12) 

(13) 


where  fi  =  blr  and  b  is  the  radius  of  the  area  covered  by  PTFE  around 
the  contacting  fibers,  see  Fig.  5.  The  total  thermal  resistance  at  each 
cell  interface  RJunc  can  be  expressed  as  a  parallel  combination  of  the 
thermal  resistances  of  2n  contact  regions. 


^junc  — 


(TCR-1  +R, 


i-1  \  1 

PTFE ' 


Fig.  4.  In-plane  SEM  image  of  a  Toray  carbon  paper,  x800  magnification. 


2  n 


(14) 
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Fig.  7.  SEM  image  of  Toray  carbon  paper  TGP-H-120  with  5%  PTFE  treatment,  x40 
magnification. 


Referring  to  Fig.  6,  the  total  resistance  can  be  written  as: 


J?tot  =  mRceii  +  (rn  -  1  )Rjunc  (15) 

where  Rcen  is  the  cell  thermal  resistance  which  can  be  found  using 
Eq.  (8). 


Rcel1  (psksWt 


(16) 


Table  2 

Input  data  for  the  in-plane  thermal  conductivity  modeling  of  Toray  carbon  papers. 


r  (p-m) 

t  (|jim) 

ks  (Wirr1  K"1) 

/+tfe  (Wirr1  K"1) 

4.25 

325  [25] 

120  [26] 

0.649  [27] 

Toray  carbon  papers  have  the  highest  through-plane  thermal 
conductivity  among  different  available  carbon  papers  with  similar 
porosity  due  to  the  contribution  of  the  binder  to  heat  transfer  in  the 
carbon  paper  GDL  [16].  The  binder  fills  the  gaps  between  fibers  and 
provides  a  better  contact.  The  thermal  conductivity  of  the  binder 
can  be  assumed  to  be  equal  to  that  of  the  carbon  fibers  [16].  The 
actual  amount  of  binders  at  contact  points  and  as  a  result  the  con¬ 
tact  area  between  contacting  fibers  are  unknown.  To  determine  the 
in-plane  thermal  conductivity,  the  value  for  a  was  estimated  as  0.1 
based  on  SEM  observation  shown  in  Fig.  4. 

The  PTFE  conductivity  is  very  low  compared  to  the  thermal  con¬ 
ductivity  of  carbon  fibers  and  its  effect  on  the  in-plane  thermal 
conductivity  is  small  as  shown  later.  Flowever,  to  include  the  effect 
of  PTFE  variation  in  the  model,  we  assume  that  the  radius  ratio  of 
PTFE  to  carbon  fiber,  /3,  is  0.25  for  5%  and  1  for  30%  PTFE  content, 
and  remains  constant  at  1  for  higher  values  of  PTFE  content.  The  fol¬ 
lowing  relationship  is  developed  for  the  PTFE  content  at  the  contact 
points  of  fibers. 

^  =  0.25  +  3(1-0.05)  (21) 


Finding  the  total  thermal  resistance  from  Eq.  (15),  one  can  eval¬ 
uate  the  effective  thermal  conductivity  using  the  following  relation. 

Kmn  =  R^Wi  (17) 

For  the  sample  with  no  PTFE  coating,  Eq.  (17)  can  be  simplified 
to: 

-  - 2 

4aT.0sfcsCOS0  ^  AodfaksCOS  0 

eff-in-°  4 aL  +  7r(m  -  1  )r  4a£cos  0  +  nr 

Based  on  scanning  electron  microscopy  (SEM)  images  of  Toray 
carbon  papers,  Fig.  7,  the  fiber  angle  with  respect  to  the  heat 
flow  direction,  0,  can  vary  arbitrarily  between  -0m  and  0m,  where 
O<0m<9O.  Considering  an  arbitrary  distribution  of  fiber  angle 
without  preferential  direction,  cos  0  can  be  estimated  through: 


cos  0  = 


Y^i=\cos®i  Y^=  icos0iA0  Y^!i=\COs0iA0 


nAO 


20  n 


(19) 


where  A 0  =  [0m  -{-0m)]ln.  The  number  of  fibers  in  a  cell  is  large 
enough,  n  =  35,000,  to  convert  the  series  in  Eq.  (19)  to  an  integral. 


cos  0  = 


Eti  cos°iA0 

20  m 


cos  OdO 
J  _ 

20m 


sin0m 

0m 


(20) 


We  measured  the  fiber  angle  in  the  SEM  image,  Fig.  7,  and  our 
analysis  shows  that  the  majority  of  fiber  angles  are  between  -75 
and  75,  i.e.,  0m  =  75.  Other  specifications  of  the  Toray  carbon  paper 
required  for  the  present  model  are  listed  in  Table  2. 


where  X  is  the  weight  fraction  of  PTFE,  0.05  <  A  <  0.3  This  is  an 
approximate  relationship,  which  is  proposed  based  on  SEM  images 
of  carbon  papers  with  different  PTFE  contents.  Flowever,  due  to  the 
very  low  thermal  conductivity  of  PTFE,  variants  of  this  distribution 
do  not  have  a  significant  impact  on  the  model  predictions. 


4.  Results  and  discussion 

Measurements  were  taken  for  TGP-H-120  samples  with  differ¬ 
ent  PTFE  contents.  A  summary  of  the  experimental  results  is  shown 
in  Table  3.  The  measurements  were  performed  at  an  average  sam¬ 
ple  temperature  of  65-70  °C.  There  is  a  small  difference  between 
the  measured  values  of  heat  fluxes  in  the  upper  and  the  lower 
fluxmeters  due  to  heat  losses  and  experimental  uncertainties.  Due 
to  heat  losses  to  the  lower  wooden  block,  the  readings  from  the 
lower  heat  flux  were  used  for  thermal  resistance  calculations.  As 
shown  in  Table  3,  the  in-plane  thermal  conductivity  appears  to 
increase  slightly  with  PTFE  content  as  a  result  of  reduced  con¬ 
tact  resistance  between  fibers.  However,  this  variation  is  within 
the  uncertainty  band  of  the  present  experimental  measurement. 
The  thermal  conductivity  values  obtained  lie  between  the  values 
predicted  by  Zamel  et  al.  [16]  and  those  of  the  manufacturer  [28]. 
Zamel  et  al.  [16]  reported  an  in-plane  thermal  conductivity  value 
of  10  Wm-1  k-1  based  on  numerical  simulations  at  a  temperature 
of  68  °C;  the  value  reported  in  the  manufacturer’s  data  sheet  [28]  is 
21  and  23  Wm-1  k-1  at  the  room  temperature  and  100°C,  respec¬ 
tively. 


Table  3 

Summary  of  experimental  results  for  Toray  carbon  paper  TGP-H-120  with  different  PTFE  contents. 


PTFE  content  (%) 

i  I  (mm) 

QuP  (W) 

Qlow  (W) 

R|ayer(KW-') 

RendtKW1) 

Keffjn  (Wm-'W-1) 

11.02 

5.72 

5.29 

63.21 

14.27 

17.39 

5 

14.34 

5.49 

5.17 

77.79 

10 

11.20 

5.61 

5.35 

65.33 

15.42 

17.33 

14.29 

5.53 

5.14 

78.88 

20 

11.37 

5.73 

5.41 

67.13 

17.19 

17.58 

14.55 

5.46 

5.25 

81.41 

30 

11.28 

5.77 

5.34 

66.95 

18.05 

17.81 

14.67 

5.62 

5.20 

82.31 
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Fig.  8.  Comparison  of  the  present  experimental  data  and  analytical  model  for  in¬ 
plane  thermal  conductivity  of  Toray  carbon  paper  TGP-H-120  over  a  range  of  PTFE 
content. 


PTFE  contents.  Toray  carbon  papers  TGP-H-120  with  PTFE  con¬ 
tent  of  5-30%  were  used  in  the  experiments.  The  experiments 
were  complemented  by  a  compact  model  for  the  in-plane  thermal 
conductivity  that  accounts  for  heat  conduction  through  randomly 
oriented  fibers,  contact  area  between  fibers,  and  PTFE  covered 
regions.  The  model  predictions  are  in  good  agreement  with  exper¬ 
imental  data  over  a  range  of  PTFE  content. 

An  important  finding  is  that  the  in-plane  effective  thermal  con¬ 
ductivity  remains  almost  unchanged,  k^l7.5Wm_1  K-1,  over  a 
wide  range  of  PTFE  content;  this  value  is  approximately  12  times 
higher  than  the  through-plane  conductivity.  However,  the  ther¬ 
mal  contact  resistance  and  the  end  effects  increase  with  the  PTFE 
content  due  to  increased  number  of  PTFE  coated  fibers. 

In  addition  to  providing  for  the  first  time  through-plane  effective 
conductivity  data,  this  work  clarifies  the  effect  of  PTFE  content  on 
the  effective  thermal  conductivity  and  contact  resistance  of  GDLs, 
and  provides  input  data  for  fuel  cell  models. 
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By  increasing  PTFE  content,  the  number  of  PTFE  coated  fibers  in 
contact  with  the  grooves  increases.  This  results  in  a  higher  contact 
resistance  between  the  GDL  samples  and  the  groove  walls.  This 
resistance  Rend  increases  26%  by  increasing  PTFE  content  from  5% 
to  30%  as  shown  in  Table  3. 

The  developed  analytical  model  for  the  in-plane  effective  ther¬ 
mal  conductivity  is  compared  with  the  experimental  data  in  Fig.  8. 
Good  agreement  is  obtained  with  a  maximum  deviation  of  5%.  We 
note  that  the  model  reproduces  the  slight  increase  in  effective  ther¬ 
mal  conductivity  with  PTFE  content  observed  in  the  experimental 
data. 

In  practice,  due  to  the  electrochemical  reaction,  overall  water 
transport  and  phase  change,  water  in  both  vapor  and  liquid  form 
is  present  in  a  fuel  cell  [2].  This  water  may  impact  the  thermal 
conductivity  while  passing  through  the  GDL.  The  thermal  conduc¬ 
tivity  of  humidified  gasses  or  water  is  several  orders  of  magnitude 
lower  than  the  thermal  conductivity  of  carbon  fibers,  and  consider¬ 
ing  the  relatively  parallel  paths  for  the  heat  transfer  through  them, 
the  effect  of  water  content  on  the  in-plane  thermal  conductivity  of 
GDLs  is  expected  to  be  minimal  when  the  GDL  is  not  significantly 
flooded.  In  the  through-plane  direction,  the  effect  of  water  is  likely 
much  more  important  as  it  can  provide  additional  pathways  for 
heat  transfer  in  the  contact  regions  between  fibers.  Burheim  et  al. 
[29]  compared  the  through-plane  thermal  conductivity  of  dry  and 
humidified  GDLs  and  showed  that  the  conductivity  increases  by 
about  70%  for  low  contact  pressures  when  water  was  added  to  the 
GDL.  This  issue  needs  to  be  further  investigated. 

5.  Summary  and  conclusions 

A  new  thermal  measurement  technique  was  developed  to 
measure  the  in-plane  thermal  conductivity  of  GDLs  for  various 
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